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Abstract 

We say a unitary operator acting on a set of qubits has been compiled if it has been ex- 
pressed as a SEO (sequence of elementary operations, like CNOTs and single-qubit op- 
erations). SEO's are often represented as quantum circuits. arXiv:quant-ph/0702144| 
by Farhi-Goldstone-Gutmann has inspired a recent flurry of papers, that propose 
quantum algorithms for evaluating NAND formulas via quantum walks over tree 
graphs. These algorithms use two types of unitary evolution: oracle and non-oracle. 
Non-oracle evolutions are independent of the NAND formula input, whereas oracle 
evolutions depend on this input. In this paper we compile (i.e., give explicit SEOs 
and their associated quantum circuits for) the oracle and non-oracle evolution opera- 
tors used in some of these NAND formula evaluators. We consider here only the case 
of balanced binary NAND trees. Our compilation methods are based on the CSD 
(Cosine Sine Decomposition), a matrix decomposition from Linear Algebra. The CS 
decomposition has been used very successfully in the past to compile unstructured 
unitary matrices exactly. 
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1 Introduction and Motivation 



(The picture on the front page is an incidence matrix for a binary tree. Empty squares 
indicate zero entries and bullets indicate non-zero entries.) 

Recently, Farhi-Goldstone-Gutmann proposed in FGG07 jl] a quantum algo- 
rithm for evaluating certain boolean formulas cj) : {0, 1}^'"'* — {0, 1} on Niys variables. 
FGG07 considers only those 0() representable as balanced binary NAND trees. These 
are perfectly bifurcating trees, with a NAND at every non-leaf node, with Ni^s leaves 
and \og2{Niys) levels, where each leaf represents a different input variable of 0(). The 
FGG07 algorithm evaluates these balanced binary NAND formulas in time 0{y/Nu}s)- 
The FGG07 algorithm uses a continuous time quantum walk. The input to the NAND 
formula is entered into the quantum walk by means of an "oracle Hamiltonian" , which 
is applied constantly and continuously during the quantum walk. One describes this 
situation as a continuous querying of the oracle. 

A few days after FGG07 appeared at ArXiv, Cleve et al published Cle07 [2] 
at the same library. This short two page paper points out that one can translate 
the FGG07 algorithm to a counterpart algorithm. The counterpart algorithm enters 
the input into the quantum walk during a finite number of interruptions or queries. 
This alternative method of entering the input is described as a discrete querying 
of the oracle. Cle07 shows that the counterpart algorithm will yield a solution to a 
balanced binary NAND formula after (9((A'";„s) 2+^) queries, for any e > 0. 

Next, Chi07 [3J, by Childs et al, appeared on ArXiv. Chi07 proposes an 
FGG07-inspired, discrete-queries algorithm for evaluating arbitrary (not necessarily 
balanced or binary) NAND formulas. The Chi07 algorithm, like the Cle07 one, re- 
quires (9((Ai';t,s) 2+*^) queries to evaluate a NAND formula. It uses a quantum walk 
on a tree graph, with a tail of nodes attached to the root node of the tree. In con- 
trast, the FGG07 algorithm attaches a runway to the tree. The runway is a finite 
line of nodes attached at its middle node to the root node of the tree. The Chi07 
algorithm uses quantum phase estimation to deduce the value of the NAND formula. 
In contrast, FGG07 uses a scattering experiment to deduce this value. 

Next, Amb07|l], by Ambainis, appeared on ArXiv. Amb07 describes an 
FGG07-inspired, discrete-queries algorithm that can evaluate balanced binary NAND 
formulas with 0{{Niys)'2) queries. Thus, the Amb07 algorithm improves over Cle07 
and Chi07 by removing the e. Amb07 also gives an algorithm that allows one to 

evaluate arbitrary NAND formulas with 0{{NiysY V'°g2 i^ivs ) queries. 

We say a unitary operator acting on a set of qubits has been compiled if it has 
been expressed as a SEO (sequence of elementary operations, like CNOTs and single 
qubit operations). SEO's are often represented as quantum circuits. 

This paper, like Cle07, Chi07 and Amb07, considers FSS07-inspired algorithms 
with discrete querying. Such algorithms use two types of unitary evolution: oracle 
and non-oracle. Non-oracle evolutions are independent of the NAND formula input, 
whereas oracle evolutions depend on this input. The goal of this paper is to compile 
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these oracle and non-oracle evolutions. We present explicit SEO's and their associ- 
ated circuits, for these evolutions. Such circuits will be required for most physical 
implementations of the algorithms proposed in Cle07, Chi07 and Amb07, but are not 
given by those papers. Such missing circuits constitute a large gap in our practical 
knowledge that this paper is intended to fill. 

Some papers (for example, Refs.[5l[6l[7l[8]) have previously addressed the issue 
of compiling quantum walks over trees and other graphs. Some of these papers use 
oracles to encode the graph topology. Our approach is different. We do not use any 
oracles to encode the graph topology. 

This paper considers only balanced binary NAND trees, although I think some 
of the methods of this paper can be applied to more general trees. 

Our compilation methods are based on the CSD (Cosine Sine Decomposition) f9], 
a matrix decomposition from Linear Algebra. This decomposition is explained in Sec- 
tion 13.41 For now, suffice it to say that the CS decomposition has been used very 
successfully in the past to compile unstructured unitary matrices exactly. 



2 Notation 



In this section, we will define some notation that is used throughout this paper. For 
additional information about our notation, we recommend that the reader consult 
Ref . |T0] . Ref . [10] is a review article, written by the author of this paper, which uses 
the same notation as this paper. 

We will often use the symbol A^^ for the number (> 1) of qubits and Ns = 2^^ 
for the number of states with A^^ qubits. The quantum computing literature often 
uses n for Nb and A^ for A^^, but we will avoid this notation. We prefer to use n for 
the number operator, defined below. 

Let Bool = {0, 1}. As usual, let Z, M, C represent the set of integers (negative 
and non-negative), real numbers, and complex numbers, respectively. For integers a, 
b such that a < b, let = {a, a + 1, ... 6 — 1, 6}. For F equal to Z or M, let F^" and 
Y>o represent the set of positive and non-negative F numbers, respectively. For any 
positive integer k and any set 5", let 5**^ denote the Cartesian product of k copies of 
S; i.e., the set of all fc-tuples of elements of S. 

We will use 9{S) to represent the "truth function"; 9{S) equals 1 if statement 
S is true and if 5* is false. For example, the Kronecker delta function is defined by 
6y = 5ix,y) = eix = y). 



EiV 



Nb- 



Let 



1 and 1 = 0. If a 
Conversely, a = 



a. 



- O-Nb- 

bin{a) 



where G Bool, then dec{a) 



We define the single-qubit states |0) and |1) by 



|0) 



1 






1 



If a e Bool'^'^ , we define the A"5-qubit state \a) as the following tensor product 



3 



(2) 



For example, 



|01) 



" 1 " 




" " 




1 







1 












(3) 



When we write a matrix, and leave some of its entries blank, those blank 
entries should be interpreted as zeros. 

Ik and Ofc will represent the k x k unit and zero matrices, respectively. For any 
matrix A G C^^^, A* will stand for its complex conjugate, A^ for its transpose, and 
A"^ for its Hermitian conjugate. For any two same-sized square matrices A and B, we 
define the o-dot product © by A S = ABA^. 

For any matrix A and positive integer k, let 



^A 



A0A, 



k copies of -A 

A®'' = A ® • • • ® A © ^ . 



k copies of -A 



(4) 
(5) 



Suppose P e Zo^ATs-i and M is any 2x2 matrix. We define M{(3) by 
M{l3) = ■ ■ ■ ® h® M ® ■ ■ ■ ® h , 



(6) 



where the matrix M on the right hand side is located at qubit position j3 in the tensor 
product of Nb 2x2 matrices. The numbers that label qubit positions in the tensor 
product increase from right to left (•<—), and the rightmost qubit is taken to be at 
position 0. 

The Pauli matrices are 



(7x = 



1 

1 



(7y = 



-i 

1 



1 
-1 



(7) 



Let a — {ax, cry, (7z)- For any a e R^, let a^ — o - a. 

The one-qubit Hadamard matrix H is defined as: 



E 



1 

71 



1 1 
1 -1 



(8) 



The iVs-qubit Hadamard matrix is defined as H®^'^ . 

The number operator n for a single qubit is defined by 



n 




1 



1 - CTz 



(9) 
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Note that 



n|0) = 0|0) = , n\l) = 1|1) 
We will often use n as shorthand for 



n = 1 — n 



1 




l + az 



Define Pq and Pi by 



Po^n 



1 




= |0)(0| , Pi=n = 




1 



Two other related 2x2 matrices are 





1 



1 




= |o)(i|. 



= |1)(0| , a+ = 

Pq and Pi are orthogonal projection operators and they add to one: 

PaPb — S{a, b)Pb for a, 6 e Bool , 



For a e Bool^'', let 



P-^ P 

^ a ^ aNg-i 

For example, with 2 qubits we have 



Pa2 ® Pa, ® Poo 



Note that 



Poo^ Po^Po^diag{l,0,0,0) , 



Poi^ Pq ^ Pi ^diag{Q, 1,0,0) , 



Pio = Pi^Po = diag{0,0,l,0) , 



Pii = Pi Pi = diag{0, 0, 0, 1] 



PgPt = 6{d, b)Pr for a,be Bool^^ , 



^ Pg = /a ® /2 ® • ■ • ® /2 = /; 



2"B 



SeBool^B 



If a G Bool^^ , and /? G Zq^at^-i, then we will denote gBooI by replacing 



aieBool 02,01,00 



ai,-,ao 



the symbol in by a dot. For example, ^ 

Next we explain our circuit diagram notation. We label single qubits (or qubit 
positions) by a Greek letter or by an integer. When we use integers, the topmost qubit 
wire is 0, the next one down is 1, then 2, etc. Note that in our circuit diagrams, time 
flows from the right to the left of the diagram. Careful: Many workers in Quantum 
Computing draw their diagrams so that time flows from left to right. We eschew 
their convention because it forces one to reverse the order of the operators every time 
one wishes to convert between a circuit diagram and its algebraic equivalent in Dirac 
notation. 

h.c. will stand for Hermitian conjugate. For example, A + h.c. will denote 



A + A\ and 



h.c. 



A 



will denote 



At 



A 



This notation is useful when A is a 



long expression that we do not wish to repeat. 

SVD will stand for Singular Value Decomposition. An SVD of a matrix A 
consists of matrices [/, V, A such that A = U where the matrices [/, V are unitary 
and A is a non-negative diagonal matrix. We will say that this SVD is "one-sided" 
if one "side" equals the identity matrix; i.e., iiU = 1 or = 1. 

When we use ni=i i,2Qi non-commuting operators Qj, we mean QiQa, 
not Q1Q2Q3 or any other permutation. Thus, the multiplication must be performed 
in the order indicated. 



3 Review of Some Prerequisites 
3.1 Circulant Matrices 

In this section, we will review some properties of circulant matrices These prop- 
erties will be used later on, to compile the evolution operator for a random walk on 
a closed loop. 

A circulant matrix is any matrix of the form 



C 



Co Cl C2 
Cn-1 Co Cl 
Cn-2 C„_i Co 



Cn-2 Cn-1 
Cn— 3 Cn— 2 
Cn-4 C„_3 



(23) 



C2 C3 C4 ■ ■ ■ Co Cl 
Cl C2 C3 ■ ■ ■ C„_i Co 

for n = 1, 2, 3, . . . and cj G C. Note that each row of C is a cyclic shift of the row 
above it. 

If we denote the eigenvalues and eigenvectors of C by 



Cv = Xv , 



(24) 
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then it is easy to show (see Ref.|TTj) that the n eigenvalues and corresponding eigen- 
vectors of C are given by 



n-l 



^m^k ^(m) 



k=0 



for m G Zo^n-i- Define matrices U and D by 



1 



UJ 
m\2 



where = e * " 



(25) 



U = (t;(0), ... ^ ^(-D) ^ jj = dtag{Xo, Ai, A2, ■ " " , A„^i) . (26) 

Then 

CU = UD , C = UDU^ . (27) 

The matrix U is in fact the discrete Fourier transform (DFT) matrix. 

In this paper, we will consider the case that ci = c„_i = g and all other Ck are 
zero. In this case, 

A„ = giu"' + {u'^r-') = giu"' + u;-^) = g2 cos( ) . (28) 

n 

3.2 Trotter Rescaling, Lie and Suzuki Approximants 

In this section, we will review the Trotter rescaling of the Lie and Suzuki approximants 
of e^"^^. More information about this topic may be found in Ref. [12]. 
Suppose A,B e C"^" and t e M. Define 

Li(t)=e*V^. (29) 

In this paper, we will refer to Li(t) as the Lie first-order approximant of e*(^+^). 
We call it a first-order approximant because, for small t, according to the Baker- 
Campbell-Hausdorff expansion [T^ [T^ . 

But what if t isn't small? Even when t is not small, one can still use the Lie 
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approximant to approximate e^^^~^^\ Indeed, if is a very large integer, then 



-La _L r 



N 



jiA+B)+4r[A,B]+0{4T) 



N 



O(-) . 



(31a) 
(31b) 

(31c) 
(31d) 



Henceforth, will refer to this nice trick as a Trotter rescaling of an approximant (in 
this case, the Lie approximant). See FigJH 




A + B 



Figure 1: Lie algebra "Physicist's picture" of Trotter rescaling. The system moves 
from to A + 5, by moving in small increments in the A and B directions. 



Next define 

S^it) = e*4e*^e*4 . (32) 

We will refer to 5*2 (t) as the Suzuki second-order approximant . One can show 1 13] 
that for small t: 

Suzuki also defined higher order approximants based on S2{t). For A; = 1, 2, 3, . . . , 
define the Suzuki {2k + 2)th-order approximant S2k+2{t) by 

S2k+2it) = Slkia2kt)S2kiil - 4:a2k)t))Sif,ia2kt) , (34) 
for some a2k ^ I^- It is possible to show that for k = 1,2,3,... and small t: 

S2k+2{t) = e*(^+^) + 0{t''+') , if a2k = ^ . (35) 

4 — 42fc+i 
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{(i2k)k=i 2 3 is a monotone decreasing sequence with 02 = 0.4145... and linijt^oo Ci2k = 
1/3. 

As with the Lie approximant, is possible to do a Trotter rescahng of the Suzuki 
approximants. One finds that for k = 1,2,3 ... , large and fixed t: 



2k+l 



(36) 



Henceforth, what we have called so far in this section will be renamed A^^, 
to distinguish it from all the other A^'s used in this paper. The T in A^ stands for 
Trotter, since it represents the number of Trotter time slices (trots?). 

If we call Nexp the number of factors of the type e"^"^ or e"^^ for some r G M, then 
we can compare N^xp and the error for the Trotterized Lie and Suzuki approximant^: 





N 


error 








'-'2k y Nt ' 


2{5)''-^Nt + 1 





Let NJI^^ stand for the number of exponentials in the Trotterized Suzuki 



exp 



approximant. Define 



^2k+l 

'nIF 



(37) 



5^ 



1 

e2fc 



This expression for Nt as a function 



Eq. ( |371) for e can be inverted to get Nt 

of e can be substituted into the value for N];^^ given in the above table, to get: 



A^, 



T-S 

exp 



2(5) 



g 2k 



+ 1 . 



(3^ 



The previous expression for A"^ has a minimum for small e. For example, suppose 



exp 



10 

T-S 
exp 



-12 



For k 

\k-l 



T-S 



1 A^ 

) exp 



g 2k 



10^ For k 



2, N^-/ ~ = 103 As k 



00, 



2(5) -^t 00. It is possible to use calculus to obtain the value of k that 



mmimizes Nl-/. One find^ 



In(^) 
21n(5) 



N 



T-S 

exp 



2t 



[21n(5)ln(i)]^ 



(39) 



^ If Nexp{S2k) is the number of exponentials in S^k, then Nexp{S2) = 3, N^xpiS/C) = 11, 



Nexp{Se) = 51. In general, Nexp{S2k+2) = Nexp{S2k)5 - 4 for A: = 1,2, 



ence equation f{k + 1) = /(fc)5 - 4 with /(I) 3 yields f{k) = 2(5 



— 0(^\k-l 



T-S 

exp 



.. Solving the differ- 
1. Finally, note that 



fik)NT-iNT-l). 



^Ref.[6j does a similar minimization for a more general Suzuki approximant that approximates 
gi(Ai+A2+...+A,„) -^jjgj-g fn> 2. We will only use to = 2 in this paper. 
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3.3 Multiply Controlled NOTs 

In this section, we will review how to compile multiply controlled U{2) operators. 
More information about this topic may be found in Ref.|15j. 

Suppose U G U{2). If r and k are two different qubit positions, gate [/(r)"*^'^^ 
(or U{t)^^'^^) is called a controlled U with target r and control k. When U = 
ax, this reduces to a CNOT (controlled NOT). If r,Ki and are 3 different qubit 
positions, ^^(t)"^''^'*"''''^"'' is called a Toffoli gate with target r and controls ki, kq. 
Suppose Nk > 2 is an integer and b G Bool^'^ . Suppose r, kjvj^_i, kj^^_2, . . . , ki, kq 
are distinct qubits and k = (k^j^^-i, ■ ■ ■ , f^i, f^o)- Gate U{T)^b(''') is called a 

multiply controlled U with target r and Nk controls k. When U = ax, this 
reduces to an MCNOT (multiply controlled NOT). 

Our goal in this section is to show a reasonably efficient way of compiling any 
multiply controlled U{2) operator; i.e., to express it as a SEO consisting of CNOTs 
and single-qubit operations. 

For any U G U{2), a multiply controlled U can be expressed in terms of 
MCNOTs as follows: 



( 

( 


1 


( 


1 

1 


( 


1 






— [ 


/ — 



|0) 



— ( 
— 1 


1 ( 

• ( 


1 — 
1 — 


— 1 


• ( 


1 — 


— 1 


• ( 


1 — 


\ 


/ ^ \ 






u 









(40) 



|0) 



This requires one ancilla qubit in state |0). The ancilla is reusable. 

Remember that our goal is to show how to compile any multiply controlled 
U (2) operator. Now all we need to show is how to compile MCNOTs. MCNOTs can 
be expressed in terms of Toffoh gates as follows: 



|0) 
|0) 
|0) 



— ( 
— < 


► i 

> ' 


► 

► 




< 




> 








— 1 




• — 












— I 


• — 








N 




1 — 






— 1 






7\ ' 




• — 


— 1 














^ 


i 



(41) 



|0) 
|0) 
|0) 



For Nk > 3 controls, this requires Nk — 2 ancilla qubits in state |0). The ancillas are 
reusable. 

So now all we need to show is how to compile a Toffoli gate. The following 
identity, true for any V G U{2), is useful in this regard: 
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1 

1 


) 

1 


— \ 


/ — 



1 


1 


• 


1 


• 1 

/ 


• 










— — v-^ — — 



(42) 



When V = Gxi this becomes: 



(43) 





















ax'' 




Ox' 




Ox' 



Thus, we can compile a Toffoh gate if we can compile a singly controlled 
(cTx)^. A singly controlled ?7(2) operator also arose in Eg. (140 p . So we should show 
how to compile a general singly controlled f/(2) operator. This can be done using the 
following quaternion trick[T6]. Observe that an arbitrary V G f/(2) need not have 
unit determinant. For example, ox and [px^' don't. So we need to express [/ as a 
product of an SU{2) matrix and a phase: U = e^'^e^^"^ , where n is a unit vector in M^, 
and 0, ^ G M. Then use the fact that any special unitary matrix e^^"^ can always be 



expressed as cxacrr, where a, 6 are unit vectors in W" . Indeed, e 



ieo 



cos 6* + in sin 6*, 



and (Tacrg 



a ■ b + ia X b, so set cos 6 - 

U 



a ■ b and n sin 6 = a x b. Thus 



Let R and S be SU (2) matrices that rotate a and b into the x axis: 

aa = RaxR'' , cr^ = SaxS'' . 

Then 



(44) 



(45) 



a 



U- P 



e*^"("Va(/5)"("V^(/?)"(") 



R 



R^S 



5t 



(46a) 
(46b) 
(46c) 



Thus, we have achieved our goal of showing how to compile any multiply 
controlled U{2) operator. 

Let G be any multiply controlled U{2) gate with controls, where > 1- 
A consequence of the results of this section is that G can always be expressed as a 
SEO with 0{Nk) CNOTs. 
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3.4 CS Decomposition 



In this section, we will define the CSD (Cosine Sine Decomposition) [9| and review 
how it has been applied to quantum computing prior to this paper. 

Suppose that U is an N x N unitary matrix, where N is an even number. The 
CSD Theorem stated that one can always express U in the form 



U 



Lo 
Li 



D 



Ro 

Ri 



(47a) 



where the left and right matrices Lq, Li, Rq, Ri are y x y unitary matrices, and 



D 



(47b) 



-Doo = Du = diag{Ci, C2, . . . , C, 



(47c) 



Dqi = diag{Si,S2, . . . ,Sn) , Diq = -D. 



01 



(47d) 



For dX\ i & N_, Ci = cos 9i and Si = sin 6i for some angle 6i . Eqs. (1471) can be 
expressed more succinctly as 



[/= (Lo©i:i)e'"^®®(i?o©i? 



(4^ 



where 9 = diag{9i, 62, ■ ■ ■ ,9n] 




Figure 2: Diagrammatic representation of CS decomposition, Eq. pTl) . 



Figl2] is a diagrammatic representation of the CSD. Note that: (l)Matrix U 
is assigned to the incoming arrow. (2)Matrix D is assigned to the node. (3)Matrices 
Rq © Ri and Lq © Li are each assigned to an outgoing arrow. 

The CS decomposition was first used for quantum compiling in Tuc99[I7]. In 
the Tuc99 compiling algorithm, the CSD is used recursively. A nice way of picturing 
this recursive use of the CSD is to represent each CSD application by a node, as 

■^Actually, this is only a special case of the CSD Theorem — the case which is most relevant to 
quantum computing. The general version of the CSD Theorem does not restrict the dimension of U 
to be even, or even restrict the blocks into which U is partitioned to be of equal size. 
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(D',x,eD'„>D>D';i') 



Figure 3: A CSD binary tree. It arises from the recursive application of the CSD. 



in Fi^ The recursion connects these nodes so as to form a binary tree, as shown 
in Figini In Figl3l we start with an initial unitary matrix Uin entering the root 
node, which we define as level 0. Without loss of generality, we can assume that the 
dimension of f/j„ is 2^^ for some Nb > 1. (If initially t/m's dimension is not a power 
of 2, we replace it by a direct sum Um © Ir whose dimension is a power of two.) We 
apply the CSD method to Uin- This yields for level a D matrix -Dg^^ two unitary 
matrices L^^ and L^^^ on the left side and two unitary matrices R^^^ and R^^^ on the 
right side. Then we apply the CSD method to each of the 4 matrices lI^\l^^\r'^^ 
and i?f ■* that were produced in the previous step. Then we apply the CSD method 
to each of the 16 R and L matrices that were produced in the previous step. And so 
on. The nodes of level Nb don't have R, L arrows coming out of them since the D 
matrices for those nodes are all 1 x 1. 

Call a central matrix either (1) a single D matrix, or (2) a direct sum Di © 
-D2 © ■ ■ ■ © -Dr of D matrices, or (3) a diagonal unitary matrix. From Figl31 it is clear 
that the initial matrix Uin can be expressed as a product of central matrices, with 
each node of the tree providing one of the central matrices in the product. We can use 
this factorization of Uin into central matrices to compile Uin, if we can find a method 
for decomposing any central matrix into a SEO. Tuc99 gives such a method. 

Let's say a few more words about central matrices. For simplicity, consider 
matrices of size 2^^ ^j^h A^^ = 3. Let Ry{(j)) = e''*"'^ . k\i S x S D matrix is of the 
form 

D= Ry{(l)ab)®Pa®Pb, (49a) 

a,b(^Bool 

for some (pab ^ An 8x8 matrix which is a sum of two D matrices is of the form 

Do©A= Pa®Ry{<P'ab)®Pb. (49b) 

a^b^Bool 

for some <^'^^ G M. An 8x8 matrix which is a sum of four D matrices is of the form 
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^00 © ^01 © ^10 © ^11 = V Pa®Pb®Ry{€b). (49c) 



a,b&Bool 



for some 0"^ G M. Thus, the D, © Di and Dqo © -^oi © Dio © -Dn matrices given 
by Egs- lHU]) only differ by a qubit permutation (apart from the fact that they have 
different angles) . Note that D of Eq. (149 ap can also be written in the form 



D = exp j i I 

\ a,b£Bool 



bCTy 0Pa®Pb] (50a) 



exp U J2 (pa,bCrY{2)Pa{l)PbiO)] . (50b) 



a,b£Bool 

For A''^ > 1, Eq. (!50|) generalizes to 



D = exp iay © ^ 0g-Pg (51a) 

\ beBool'^B-i / 

= expltariNB-l) Yl - 2, • • • , 2, 1, 0) . (51b) 

\ beBool'^B-^ / 

In general, a central matrix acting on Nb qubits is a matrix of the form exhibited by 
D in Eq. llST!) . or a qubit permutation thereof, or a diagonal unitary matrix. 

In my papers that followed Tuc99, I've begun calling a D matrix of the form 
Eq. (l5Tl) a "multiplexor"]^ When I want to be more precise, I call it an Ry{2)- 
multiplexor with target qubit Nb — 1 and control qubits A''^ — 2, . . . , 2, 1, 0. The 
Ry{2) term refers to the fact that the set of operations acting on the target qubit are 
2x2 qubit rotations Ry{4>) = e*'^'^^ for some e M. More generally, one can speak 
of f/(2)-multiplexors. Henceforth in this paper, I'll continue using this multiplexor 
nomenclature, even though it's not used in Tuc99. 

Tuc99 gives identities for decomposing an arbitrary i?y(2)-multiplexor and a 
diagonal unitary matrix into a SEO with 2^^ CNOTs. FigJH shows an example of 
the SEO decomposition found in Tuc99 for an i?j^(2)-multiplexor. In FigSJ 0,1,2,3 
are the control qubits, and 4 is the target qubit. The empty square vertices represent 
Ry{2) gates. The symbol to the left of the equal sign, the one with the "half- moon" 
vertices, was invented by the authors of Ref.[IS] to represent an i?y(2)-multiplexor. 

The Tuc99 algorithm is implemented in a computer program called Qubiter, 
available at Ref. [191. 



^"multiplexor" means "multi-fold" in Latin. A special type of electronic device is also called a 
multiplexor or multiplexer. 
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I iTVi r7vi r?vi ryvi ryvi rTsi itvi [Tvi rysi iyvi itvi ryvi r/vi rTvi fTvi 

Figure 4: A possible decomposition of an i?j^(2)-multiplexor with 4 controls. 



Later papers [201 CSl EI] have improved the Tuc99 algorithm. Instead of grow- 
ing the CSD tree all the way down to where the leaf nodes contain diagonal unitary 
matrices, they stop one level earlier, when the leaf nodes contain ?7(2)-multiplexors. 
Then they apply a nice technique, due to Vidal and Dawson[22], for expressing ar- 
bitrary U{2) matrices using just 3 CNOTs. They are able to express an arbitrary 
structureless [7(2^'^) matrix using ^4^'^ + 0{2^^) CNOTs (and also some qubit rota- 
tions, of course). It's easy to show[53] that the number of CNOTs needed to express a 
structureless U{2^^) matrix is greater or equal to |(4^-s — ?>Nb — 1)- Roughly speak- 
ing, in these algorithms, a unitary matrix is decomposed into 2^-s (2) -multiplexors, 
and each of these multiplexors is expressed as a SEO with 2^^ CNOTs. 

Note that quantum compiling papers usually measure circuit complexity by 
counting the total number of CNOTs. Others count the total number of elementary 
gates. One could also count just the total number of control vertices in the circuit 
diagram. These 3 measures are linearly related so they have the same big O behavior 
when expressed as a function of the number of qubits. It is expected that the quantum 
computers of the future will be able to perform single-qubit operations much faster 
than CNOTs, so it makes sense to ignore the single-qubit operations and count only 
the CNOTs. 

The Tuc99 algorithm and its variants expect a structureless unitary matrix as 
input. They do not assume any symmetries in the structure of the input matrix being 
decomposed. But the CSD does have some free degrees of freedom (for example, 
certain phase choices) which are chosen arbitrarily. The number of these degrees 
of freedom increases as the symmetry of the input matrix increases. If the input 
matrix does have a symmetric structure, it is sometimes possible to chose those free 
degrees of freedom in a way that reduces the length of the output SEO. For example, 
Refs. |24l [25] show that the Coppersmith decomposition ^26] of the discrete Fourier 
transform matrix can be obtained using the Tuc99 algorithm and some judicious 
choices for the free degrees of freedom. 
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4 Some Identities for 

Exponentials of Partitioned Matrices 

In this section, we present some identities satisfied by exp(iM), where M is a Hermi- 
tian matrix partitioned into four equal sized blocks. 

Suppose Moo, Mqi, Miq and Mu are square matrices of the same dimension, 

and 



M 



(52) 



Moo Moi 
Mio Mu 

We will say M is block diagonal if Mio = Moi = and block anti-diagonal if 

Mqo = Mil = 0. We are interested in the case were M is Hermitian. 

If the M of Eq.([52]) is block diagonal, then exp(iM) = exp(iMoo) ©exp(iMii). 
But suppose M is block anti-diagonal. What can we say about exp{iM) then? 

Lemma 1 Suppose F G C"'^'^ and F = VAW is an SVD, so V, U are unitary 
matrices and A is a non-negative diagonal matrix. Then 



exp(z 



" 


Ft " 




" u 


" 




F 





) = 





V 





cos(A) isin(A) 
zsin(A) cos(A) 



W 





(53) 



proof: Let LHS and RHS stand for the left and right hand sides of Eq. (l53!) . Then 

(54a) 



LHS 



exp[t 

U 
V 

RHS . 



' u 


" 




' 


A " 




■ f/t 








V 




A 










_ 



exp[t 



A 
A 



f/t 

yt 



(54b) 
(54c) 



LHS can still be exponentiated without having to find the SVD of F. See 
Appendix |X] if interested. 
QED 

Suppose the M of Eq. (l52l) only has Mn = 0, so it isn't fully block diagonal. 
What can we say about exp(iM) then? 



Lemma 2 Suppose A, B ^ 



and A is a Hermitian matrix. Then 



exp(z 



A 5t 
B 



exp(z 



A 




|exp(2^ 



dt 



h.c. 

Q^itA Q 



(55) 



where T{} indicates a time-ordered exponential. 
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first proof: Define the following three matrices: 

Q 

For any t G M, consider the t-dependent matrix 

x{t) 



' A 


5t " 


! Qa — 


'AO' 


, Qb = 


' 


5t " 


B 








B 






The derivative of x{t) is 



(56) 



(57) 



x'{t) = e-'"^^{-iQA + tQ)e'"^ 



Hadamard's identity|l3j tells us that 



e-**^-^ Qb e 
It is easy to check that for n = 1, 2, 

[Qa, TQb 

Therefore, 



HQa — p-it[QA, ] 



Qb. 



A"5t 
fi(-A)" 



(58a) 
(58b) 



(59) 



(60) 



ie~'^^^ Qb e'^^^ 



E 



-ity 



n=0 

e-^*^5t 

^^itA Q 



v4"5t 
B{-AY 



Thus, x{t) defined by Eq. (l57|l satisfies the differential equation: 



e-^*^5t 

^^itA Q 



x'{t) = i 

This differential equation for x{t) is of the form 

x'it) = M{t)x{t) 



x{t) 



(61a) 
(61b) 



(62) 



(63) 



where M{t) is a t-dependent matrix. If we express the derivative on the left hand 
side of Eq. (l63|) in terms of infinitesimals, we find 



(64) 



x{t + dt) ^ [1 + M{t)dt\x{t) ^ e*^(*)'^*x(t) 
Iterating Eq.l lMll gives: 
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= T|e^*^*'*^(*')}x(0) . (65) 

Eq. flBSl) is the solution of any differential equation of the form Eq. flB^ . We can use 
Eq. (!65l) to solve Eq. (l62!) with the initial condition x(0) = 1 (from Eq. (!57|l ). 

second proof: Here is an alternative proof that uses infinitesimals instead of 
differential equations. For any large integer N, let At = and tk = kAt for k G ^o,Af- 
Thus, to = and t^r = 1. Then 



^ ^-iQA/Qi^tQB^iAtQA\N 



n {< 



-itjQ A J.AtQ B „itjQ A 



} 



i=Af,..,2,l 

W exp(e-^*^«-^ iAt Qb e**^^-^) 



j=Af,..,2,l 



JJ^ exp(iAt 

i=Af,..,2,l 

T < exp(i / dt 



e-**^^5t 



e 



-itA 



5t 



(66a) 
(66b) 

(66c) 
(66d) 
(66e) 



QED 

In the previous Lemma, the magnitude of the matrices A and B was arbitrary. 
If we assume that A and B are 0{g), where g is small, then one can say more about 
the exponential of an M such that Mn = 0. 

Lemma 3 Suppose A,Be C"^", A is a Hermitian matrix, and A,B are both 0{g), 
where g is small. Then 



exp(z 

where 



A fit 
B 



) = exp(i 



A 




exp( 



Ti + T2 iB^ 



iB 



: BAB'f 



+ 0{g') , (67) 



Ti = —{B^BA + AB^B) , 
T2 = {^){-B^BA^ + A^B^B) 

— A -A 

B = B sinc(— )e*2 . 



(68) 
(69) 
(70) 



Also, to lower order, 



18 





" A 


5t " 




" A 


" 









exp(i 


B 





) = exp(z 








) exp(z 








0{g') . (71) 



proof: (The arXiv source code for this paper includes an Octave/Matlab subroutine 
called qtree_cbh.m that checks this Lemma.) If x and y are both of order g and g is 
small, then the Campbell-Baker-Hausdorff Expansion is [131 



11 1 

e"e^ = exp <{ X + y + -[x, + — ([x, .fy + [y, .fx) + — [y, [x, [y, x]]] + 0{g''] 



(72) 



If and Qb are defined by Eq. fl56|) . then this Lemma follows from Eq. f[72l) with 



X = -iQa, y = KQa + Qb) ■ (73) 
The definition of B was inspired from Lemma [2] and the observation that 

"1 



itA 



dt e 



where sinc(x) = sin(x)/x as usual. 
QED 



smcl — le 2 
iA ^2' 



(74) 



5 General Strategy for 

Compiling Any Quantum Walk 

In this section, we describe our general strategy for compiling the FGG07 algorithm 
(or, more precisely, the FG GO 7- inspired algorithm with discrete queries that was first 
proposed in Cle07). The compilation strategy described in this section is also useful 
for compiling other quantum walks. 

The standard definition of the evolution operator in Quantum Mechanics is 
U = e~**^, where t is time and if is a Hamiltonian. Throughout this paper, we will 
set t = — 1 so f/ = e*^. If H is proportional to a coupling constant g, reference to 
time can be restored easily by replacing the symbol g by —tg, and the symbol H by 
-tH. 

FigOl^a) shows an example of a "loop" graph. It has 8 nodes (nodes=states), 
labelled 0, 1, . . . , 7. We wish to consider the following Hamiltonian for transitions 
along the edges of this loop graph: 
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(a) 



(b) 



5 6 

(c) 



Figure 5: (a)loop graph (b)tree graph (c)tree and loop graphs glued together 



Hi 



ip 








1 


2 


3 


4 


5 


6 


7 







1 












1 


1 


1 




1 












2 




1 




1 










3 






1 




1 








4 








1 




1 
















1 




1 




6 












1 




1 


7 


1 












1 





(75) 



where G M. 

A binary tree with A + 1 levels has 1 + 2 + 2^ + . . . 2"^ = 2^+^ - 1 nodes. To 
reach 2^+^ nodes, we add an extra "dead" or "dud" node, labelled with the letter 
"(i". This d node is not connected to any other node in the graph. If we include 
this dud node, then the number of leaves Ni^g is exactly half the number of nodes: 
Ni.s = |2^+i = 2^. 

FiglSl^b) shows an example of a binary-tree graph with A = 2. It has 8 nodes, 
labelled d, 1,2, ... ,7. We wish to consider the following Hamiltonian for transitions 
along the edges of this binary-tree graph: 
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(76) 



where gp^q G M for all p, q. 

Yig^c) shows an example of the graph used in FGG07. It contains a binary 
tree attached to a loop which serves as a "runway" for a wave packet. For the 
algorithm of FGG07 to work, the runway must be much longer than ^/Ni^s as A^^^ — > 
oo|j In Fig|5](c), the tree has 8 nodes (including the dud one), and the runway has 8 
nodes. Let's call the tree nodes tr^, where j G Zqj, and the runway nodes Ipj, where 
j G Zqj. We wish to consider the following Hamiltonian for transitions along the 
edges of this tree and runway graph: 



Hip 


























hgl 








Htr 












"'in 

















(77) 



If g is the coupling constant that appears in the loop Hamiltonian Hip^ then we will 
assume, as FGG07 does, that the coupling constants {gp,q\'ip,q that appear in the 
tree Hamiltonian Htr are all equal to g. Besides Htr and iJ;p, which have already 
been discussed, the Hamiltonian of Eq. (!77|) includes blocks with hgi, hin, and their 
Hermitian conjugates. The tree-loop "glue" Hamiltonian Hgi is defined by 

Hgi = hgi + = g{\lpo){tri\ + h.c.) . (78) 

It corresponds to the graph edge that connects the tree and the runway. The input 
(oracle) Hamiltonian Hin is defined by 

Hin = hin + = g^Xk{\trk){ink\ + h.c.) , (79) 

k 

where k runs over all tree leaves, and Xk G Bool are the NAND formula inputs. 

^ FGG07 "requests" a runway of length Ni^s- However, only a stretch of length c^/ Ni^s , for 
some large constant c, is ever used(traversed) by the narrow wavepacket during the duration of the 
experiment. Chi07 and Amb07 request a tail of length c\/Nivs- 
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The compilation strategy of this paper is to spht the full Hamiltonian of 
Eq. (1771) into two parts, which we call the bulk Hamiltonian H^uik and the bound- 
ary corrections Hamiltonian H^orr'- 

H = Hbulk + Hcorr ■ (80) 

We will set 

Hbulk = Hip + Htr , Hcorr = Hgi + • (81) 

Note that [Hip, Htr] = and [Hgi, Hin] = 0. In subsequent sections we will show how 
to compile the evolution operators exp{iHip) and exp{iHtr), and therefore exp{iHbuik)- 
We will also show how to compile exp{iHgi) and exp(zifj„), and therefore exp{iHcorr)- 
Exact compilation of exp{iHbuik + iHcorr), although possible via an algorithm such as 
that of Tuc99, is numerically very laborious. Our strategy is to avoid such numerical 
calculation by combining the compilations of exp{iHbuik) and exp[iHcorr) via the 
Trotterized Suzuki approximation. 

Actually, FGG07 assumes a runway that is a straight line, but replacing a 
straight line by a large enough closed loop does not change the algorithm significantly. 
Purists can always convert the loop Hamiltonian into a line Hamiltonian by adding 
an extra loop-cutting Hamiltonian 

Hip-cut = -g{\lpz){lp4\ + h.c.) (82) 

to the corrections Hamiltonian. It is also possible to convert a loop into a "tail" like 
the one used in Chi07 and Amb07, by adding a loop-cutting Hamiltonian that cuts 
the loop at the edge that connects nodes IpQ and Ipj. 

6 Compiling Loop Graphs 

In this section, we will show how to compile (exactly) the evolution operator e^^^^ for 
a loop graph. 

Eq. 0751) is an example of a loop Hamiltonian Hip. Since Hip is a circulant 
matrix, e^^'-p = Ue~'^^W , where U is a. discrete Fourier transform matrix and D is a 
diagonal real matrix. The U and W can be immediately compiled via the Copper- 
smith decomposition|26], using O(A^B^) CNOTs, where A^^ = 2^^ jg the number of 
states in the loop. e~*^ is a diagonal unitary matrix. Tuc99 shows how to compile 
an Ns X Ns diagonal unitary matrix using 0{Ns) CNOTs. 

Suppose 9j; eM. for all b G Bool^'^ . Suppose (3 are Nb distinct qubit positions 
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and a is the position of an additional ancilla qubit. Note that 



exp [lY^hPifi) I = exp I 2az(«)X]%Pg(/5)j |0), (83a) 

-if<^x(a)gxp I 2ay(a) J]%Pg(/3) e^^"^(")|0)„ (83b) 



e 



We used o"^|0) = |0) and e^^^"^ OYe"^"^ = ctz- Thus, if one can compile any Ry{2)- 
multiplexor with Nb controls, then one can compile any 2^^ x 2^^ diagonal unitary 
matrix. 



7 Compiling Glue and Cuts 

In this section, we will show how to compile the evolution operators e*^™' for a loop 
cut and e*^f' for the tree-loop glue. 

Consider the loop cut first. For example, suppose we want to cut the loop 
graph of Figj5]^a) at the edge connecting the states |0) and |7). 

Define 

U=[ax{2)ax{l)r^'K (84) 

Then 



\7><0\ + h.c. = |lll)(000| + /i.c. (85a) 
= ?7(|001)(000| + /i.c.)[/"^ (85b) 
= t/(Tx(0)^(2)«(i)f/t . (85c) 

To get Eq. fl85cl) . we used the fact that <Jx{P) = (|0)(1| + h.c.)i3 for any qubit position 
j3. Finally, note that 

In general, suppose we want to cut the edge that connects states \j) and \k), 
where j and k are the decimal names of the states. Let j = bin{j) and let jp be the 
component of j at bit position p. Define k and kf^ analogously, j k so j and k must 
differ at one bit position, at least. Let (3o be one such bit position. Assume jp^ = 1 
and = 0. Define 

U=[ n ^xif3)r^^"^- (87) 

Then 
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\j >< k\ + h.c. = Uaxi/3o)^^'^^^''^^''^''^^^U^ ■ (88) 

Therefore, 

^iHcut _ ^-ig{\j)(k\+h.c.) _ jj]^^-ig(Tx{l3o)^f3:f3^i3^Pkfi(P) jj] _ ^gg^ 

Compiling e^^f' is identical to compiling e*^™* . The only difference is that the 
coupling constant g is replaced by —g. 



8 Compiling Line Graphs 

In this section, we will show how to compile (approximately) the evolution operator 
6*^""*= of a line graph (open loop). Previously, we showed how to compile the evolution 
operators e*^'p and e*^™* for a loop graph and a cut. A possible compilation of 6*^""*= 
can be achieved by combining e*^'^ and 6*^^="* via Trotterized Suzuki. Alternatively, 
one can compile 6*^""*= directly, as will be shown in this section. 

Eq. (|75i) is an example of a line graph Hamiltonian Hune, provided we set 
to zero the entries (0,7) and (7,0). In Eq. (|75l) . we've label states by the decimal 
numbers to 7, and we have assumed that transitions between these states are such 
that their label x G Zq j can only vary by Aa; = ±1. As in spectroscopy, let's call 
this constraint on Ax a selection rule. Although compiling e*^""^ with an Hune that 
satisfies Ax = ±1 is possible, a cleaner, simpler compilation [H] can be achieved if 
Hiine satisfies a different selection rule; namely, the constraint that states can only 
change to other states iff the initial and final states are adjacent in a Gray ordering. 
In a Gray ordering, states are labelled by a binary number, and adjacent states have 
labels that differ only at one bit position. 
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Figure 6: Line graph for eight states. States are connected iff their binary represen- 
tations differ only at one bit position; i.e., connected states are adjacent in a Gray 
ordering. 



Figl6] shows eight states, ordered linearly in a decimal ordering from left to 
right. Edges connecting the states indicate possible transitions between the states. 
These transitions can only occur between states that are adjacent in a Gray ordering. 
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The Hamiltonian Hum that describes transitions along the edges of the graph of Figj6] 
is the following. (Note that the states along the rows and columns of Hune are ordered 
in a decimal ordering.) 
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g{I ^ I ^ ax + I ® (Tx ® Pi + cTx ^ Pi ^ Po} 
g{axiO) + ax(l)n(0) + ax(2)n(l)n(0)} . 



where G M. Define 



A = gaxiO) 

B, = gaxilMO) 

B2 = gax{2Ml)n{0) 



(90) 
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(91) 



(92) 

(93) 
(94) 



(95) 



Note that the Bj commute because they live in different states (see Figj6]). Thus, 

(97) 



Rx 



+ 0{g' 



if we use a first-order Lie approximant, and 
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if we use a second-order Suzuki approximant. 
With four bits, one gets 



(98) 
(99) 
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(100) 



where g e 



= g{ax{Q) + (7x(l)n(0) + (7x(2)n(l)ri(0) + ax(3)n(2)ri(l)ri(0)} . (101) 
Define 



B, = gaxilMO) 
B2 = gax{2)n{l)n{0) 
= gax{3)n{2)n{l)n{0) 



(102) 



Then 



— e 



i{A+Bi+B2+B3) _ JAJBiiB2jB3 



— Rx 



Rx 



Rx 



+ 0{g' 



(103) 
(104) 



Rx — 



if we use a first-order Lie approximant, and similarly if we use a second-order Suzuki 
approximant. 

Generalization to arbitrary number of bits is obvious. 
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9 Generalities about Compiling Tree Graphs 



In this section, we will introduce several key ideas that are useful in compiling the 
evolution operator ex]i{iHtr) for a general binary tree graph. 

For A = 1,2,3,..., define the 2^ x 2^ matrix to be, up to a constant 
factor, the Hamiltonian Htr for a binary tree with 2^ nodes. By studying Eq. (l76p . 
one realizes that the family of matrices {yl2A}vA can be specified in terms of a family 
of matrices {B2\}\f\ as follows. Let A2 = 0, and 



A2X+1 



A2X Bl 
B2X 



(105) 



for A = 1,2,3, ■■ ■. (The subscripts of A2X and B2X indicate the dimension of these 
matrices.) For example, according to Eg. (11051) . Aiq has the form: 
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(106) 



9.1 CSD-ready Evolutions 

By applying Eq. (!7T!) repeatedly, we can factor e'"^^^ with Aiq given by Eq. (ll06p as 
follows: 



where 



(107) 



and 



F3 = exp(z 
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108a) 



;i08b) 



:i08c) 
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is defined for A G Zi 3 by 

= B2X smc(^)e^^ . (109) 

The above example generalizes as follows. For A = 1,2,3, . . ., and A G ^i,a, 
we hav^ 

e'^^^+' =T,r2...Tj, + 0{g^) , (110) 



(111) 



and 

= B2X smc{^)e'— . (112) 

Since these are ready for an application of the CSD, we will henceforth refer 
to them as CSD-ready evolutions (or factors) of exp(iy42A+i), and to Eg. (11 101) as 
a factorization of exp{iA2A+i) into CSD-ready factors. 



r> 



exp[% 
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02A+i_2^+i 



9.2 SVD of Bo 



Next we will calculate an SVD for each of the matrices B2X that appear in Eq. (llOSp . 
These SVD's will be useful to us in future sections. 

We will use E4 to denote the 4x4 matrix that exchanges (swaps) 2 qubits: 
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Given a,-, bj G M for some index j, define 



-aj bj 



The matrix B2 that appears in Eq. fllOSp is of the form 

„ a 
b 



(113) 



(114) 



(115) 



^ An important special case is when the d\ (defined by Eq. (|13ip ) are all equal to V^g. In this 
case, a slight redefinition of F"^ in Eq. (|llH) makes Eq. (|110p good to order 0{g'^) instead of 0{g^). 
See Appendix IbI 
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These values of a and b define a matrix U and a scalar p via Eg. (11141) (with the j 
subscript absent). The SVD of B2 is one-sided and can be found in closed form. It's 
given by: 




Ai 



where 



and 



Fo = U 



(116) 



(117a) 



Ai = p . (117b) 
When expressed in operator rather than matrix notation, Eqs. (lll7p become 



i^2(0) = f/(0) 



;il8a) 



and 



Ai = p. 

The matrix B4 that appears in Eq. fllOSp is of the form 



;il8b) 



Ba 



ao 
bo 



ai 
61 



(119) 



For j = 0, 1, these values of aj and bj define a matrix Uj and a scalar pj via Eq. (lll4p . 
The SVD of B^ is one-sided and can be found in closed form. It's given by: 



B. 



f/n 





Po 




Pi 



where 



and 



A. 



F4 = (f/o © Ui)E, 



(120a) 
(120b) 

(121a) 



A2 = diag{po,pi) 



(121b) 
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When expressed in operator rather than matrix notation, Eqs. fll2ip become 



F4(l,0) = f/o"(^)(0)f/r^^) (0)^4(1,0) , (122a) 

and 

A2(0) = J2 PbPbiO) . (122b) 

beBool 

The matrix that appears in Eq. OlOSp is of the form 







aoo 
boo 










aoi 
feci 














aio 
bio 










an 
6n . 



(123) 



For j G {00, 01, 10, 11}, these values of aj and bj define a matrix Uj and a scalar pj 
via Eq. flll4p . The SVD of Bg is one-sided and can be found in closed form. It's given 
by: 



00 



B. 



01 



10 



11 









Poo 












Poi 
















Pio 











pii . 



(124a) 



where 



and 



A4 



Fs = {Uoo © f/01 © Uio © f/n) (/ © ^4) (^4 © /) 



(124b) 



(125a) 



A4 = (im5((poo,Poi,Pio,Pii) • 



(125b) 
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When expressed in operator rather than matrix notation, Eqs. fll25p become 



F8(2,l,0) = f/oo(0)"(2)"«f/oi(0)"(2Wi)f/,o(o)"(2Wi)f/^^(0)«{2Wi)^^ ^ 

(126a) 

and 

A4(l,0)= Yl P?^6(1'0)- (126b) 

bf^Bool^ 

The above example generahzes as follows. For A = 1, 2, 3, ... , 

F2.(A- 1, . . . , 2, 1,0)= n {^/^'^'"''■■■''''\0)}E4(0, 1)^4(1, 2) . . . E4(A-2, A- 1) , 

(127) 

and 

A2A-i(A-2,... ,2,1,0)= Pg^g(A-2,... ,2,1,0). (128) 

In this paper, we are interested mainly in the case of balanced binary NAND 
trees. For such trees, and are both independent of b. In particular, 

= HcTx = e'^"^ (129) 

for all b. Therefore, for balanced binary trees, is proportional to the identity 
matrix and 

F2a(A-1,...,2,1,0) = (iJax)(0)E4(0,l)^4(l,2)...E4(A-2,A-l) , (130) 
for all A = 1,2,3.... 

10 Compiling Balanced Binary NAND Trees 

In Section [9l we found an SVD for each matrix i?2A. Half of its singular values 
were zero and the other half were stored in a diagonal matrix called A2A-1. In this 
section, which contains several subsections, we will give an approximate compilation 
of exp{iHtr) for a special type of binary tree that includes balanced binary NAND 
trees. The special type of binary trees that we will consider in this section satisfies 

A2A = rfA/2A , (131) 

for all possible A, where dx G M. In other words, the A2A of Htr are all proportional 
to the identity matrix. Balance binary NAND trees satisfy Eq. fll3ip because their 
coupling constants gp^g are the same for all p, q. 



1 1 
-1 1 
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10.1 Diagonalizing the Tree Levels 

For A = 1,2,3, define 



B2A = 



O2A-1 



(132) 



A2A-1 



A family of matrices {v4.2a}va can be defined in terms of tlie family of matrices {i32A}vA 
as follows. Let ^2 = 0, and 



for A = 1,2,3, The subscripts of ^2^ and ;B2a indicate the dimension of these 

matrices. 

Next, we will explain how the matrices A2X are related to their namesakes A2>^- 
For illustrative purposes, consider Aiq and Aie first. We claim that if we define 

J^i6 = diag{Is,Fs)diag{l4,F4,l4,Fi)diag{l2,F2,l2,F2,l2,F2,l2,F2) 



A2X Bl 
B2X 



(133) 



{h®F^){h®F4f\l2®F2) 




(134) 



then 



A 



16 — ^ I&AiqT iQ . 



(135) 
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This claim follows from the following observations: 



A 



16 



02 


Bl 


B\ 






O2 


Bl 


^4 


O4 




Os 



136a) 





h 


F^ 


h 


F^ 


h 


F^ 


h 




h 


Fl 


h Fl 


h Fl 


h 


h 


h 















F2 









A2 


F4 


h 










A4 


F2 




A2 




h 


h 








F2 


F, 


h 




A4 


F2 



^1 



16 



02 


Bl 


B\ 




B2 


O2 


B^ 
"s 




O4 


Bs 


Os 



•/ 1 



16 



(136b) 



(136c) 



= J-i6^i6^iV (136d) 

Going from Eq.f ll36al) to Eq. (1136bl) is the crucial step that requires assumption 
Eq. fll3ip . Eq. fll36bl) is a non-conventional way of expressing the product of 3 ma- 
trices: the matrix that is totally enclosed by double lines is being pre-multiplied by 
matrix jF^g and post- multiplied by matrix ^|g. jF^g is spread out over the left margin 
and J-'iQ over the top margin of the matrix enclosed by the double-line border. Writ- 
ing Eq. fll36bl) in this nonconventional way makes it clear why assumption Eq. (11311) is 
needed. For example, because A4 acts like a scalar, the matrix F^ in the last column 
of the top margin commutes with A4 and can cancel the F2 in the fourth row of the 
left margin. 
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Note that Tvo given by Eg. (11341) can be expressed in operator rather than 
matrix notation as follows: 



^16(3, 2, 1, 0) = F8(2, 1, 0)"(3)F4(1, 0)"(2)F2(0)"« 
The above example generalizes as follows. For A = 1,2,3..., 



(137) 



^2A+i (A, . . . , 2, 1, 0) = (A - 1, . . . , 2, 1, O)'^(^) . . . F4(l, 0)"(2)F2(0)"« , (138) 



and 



— J- 2A^2^-^2^ • 



Thus, for A = 1,2,3, 



eXp(iy42A+l) = JF2A+1 exp(M2A+i)-^2A+i 



(139) 



(140) 



10.2 Separating into Strands 

A consequence of diagonalizing the levels of the binary tree graph is that the tree 
graph is replaced by a collection of subgraphs shaped like strands. In this section, we 
will factor exp(z^2A+i) iiito separate contributions from each of these strands. 

Let us consider an example first, before dealing with the general case. Suppose 
A = 5 and consider exp(z^64)- By studying carefully the definition Eg. (11331) of A%a^, 
one can see that 



15 



exp(i^64) = ]^exp(i^i 



64 ) ) 



(141) 



s=0 



where the matrices J^^\ are defined in terms of matrices and B^^l as follows. The 
index s labels the 16 "strands" that contribute to exp(i^64)- Suppose s E ^o,i5- Let 



B. 



O2A- 



A 



(142) 



for A G Zi 5. Let 



is) 



and 



•^2^+1 



As) K?{s)t 



B. 



is) 



(143) 



for A G Z'l 5. 



For s = 0, one finds 
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•^32 

•^16 

(0) 
2 



A. 



(cii|l)(3| +d2|3)(7| + rf3|7)(15| +rf4|15)(31|) + /i.c. 
(di|l)(3|+ci2|3)(7| + rf3|7)(15|) + /i.c. 
{d^\l){3\+d2\3){7\) + h.c. 
di\l){3\ + h.c. 




For s = 1, one finds 



r(0) 
^32 

(0) 



^4 

i3f 



(0) 



4|31)(31| 
d4|15)(15| 
d3\7){7\ 
t^2|3)(3| 
d,\l){l\ 



(144) 



•^32 

A^'^ = 



(d2|2)(6| + o?3|6)(14| + d4|14)(30|) + h.c. 
{d2\2){6\+d3\6){U\) + h.c. 
{d2\2){6\) + h.c. 




*-*32 

r(1) 



4|30)(30| 

rf4|14)(14| 

d3|6>(6| 

rf2|2)(2| 





(145) 



For k e Zo,i (and s — 2 + k E ^2,3), one finds 



A- 



(2+fc) 



32 _ (rf3|5_A;)(13-A;|+4|13-/c)(29-A;|) + /i.c. B32 
{d3\5-k){13-k\) + h.c. 



•^16 



A^'"-'^ = 



^(2+fe) 



/I (2+fc) 



For k e Zo,3 (and s = 4 + /c e ^4,7), one finds 



(2+fc) 



c^5|29- A;)(29- A;| 



Bf/'^) = 4|13-A;)(13-A;| 



(2+fc) 
8 

(2+fc) 



d3\5-k){5-k\ 

^(2+fc) ^ 

(146) 



j(4+fc) 
•^32 

^(4+fc) 
•^16 

^(4+fc) 



(i4|ll-A;)(27-A;|+/i.c. H 







^(4+fe) 



(4+A:) 
32 
(4+fc) 
16 

(4+fc) 



4|27- A;)(27-A;| 
di\ll-k){ll-k\ 



^(4+fc) 



^(4+fc) 







(147) 



For k e Zqj (and s = 8 + A; e ^8,15); one finds 



/((S+fc) 
•^16 



32 

(8+fc) 





.(8+fe) 

^8 



.(8+fc) 



.(8+fcO 
^2 



w(8+fc) 
*-'32 

Bf/^) 



rf5|23 - A:) (23 - k\ 
^{8+fc) ^ ^(8+'^) = B^^+''^ 



The above example generalizes as follows. For A = 1,2,3, 





(148) 



exp(i^2A+i) = Yi exp(i^^lVi) > 



(149) 



s=0 
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where ^f,A+i is defined in terms of matrices and B'^x as follows. Suppose s G 



V2A+ 
Zq 2*^1—1- Let 



2^ 



O2A- 



(150) 



for A e Zi,A. Let A^"'^ = and 







(151) 



for A G 



3 7_ 11 15 19 23. 27 31 



11 



15 



19 



23 



27 



31 



Figure 7: Illustration of strands as entries of matrix Az2- An empty cell represents a 
zero matrix entry. The number inside a non-empty cell is the number of the strand 
that owns that cell. 



Figures El El and [7] are three alternative ways of visualizing strands. 

It is convenient to define infinitely many strands, and to define each strand 
as an infinite sequence of integers. The integers in a strand label some nodes of an 
infinite binary tree. When we consider a finite tree, the strands are truncated to a 
finite length. 
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16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 

7 6 5 4 3 2 1 

strand number 

Figure 8: A binary tree with 32 nodes, including a dud node d. Heavy blue lines 
mark the node strands of the tree. 

For s = 0,1,2, let 

A™„(s) = \\og^{s + 1)1 + 1 . (152) 

Thus, 



s 







1 


1 


2 


2,3 


3 


4-7 


4 


8-15 


5 



(153) 



By the strand s, we mean the infinite sequence 

Note from our A = 5 example that, in the definition of ^2^^, strand s is truncated 
when it reaches the term: 

z//(s,A) = 2^-5-1 . (155) 
Thus, for A G 2'i a and s G Zo 2a-i_i, 

A-l 

A;l= dy\2'' -s-l){2''+'-s-l\, (156) 

and 
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48 
















56 














63 



Figure 9: Each column of this "log2 staircase" gives a different strand of the binary 
tree in FigJHl Strand numbers are given at the top edge. 



B)^> =dx\2^ - s-l){2^ - s-l\ . (157) 

Note that 2^^^ strands ( 2^^^ = a quarter of the number of nodes and half the 
number of leaves of the tree) are active (come into play) in the definition of ^2^+1 • 
For example, 16 strands are active in the definition of ^64? and Aqa is associated with 
a tree with 64 nodes and 32 leaves. 



10.3 Factoring into CSD-ready Evolutions 

In this section, we will factor e 2A+1 for each s into CSD-ready evolutions. 

Let us consider an example first, before dealing with the general case. Suppose 
A = 5 and consider exp(2^g4'*). By applying Eq. (17T]) repeatedly, we can factor e*-^^* 
into CSD-ready evolutions, as follows: 

e^4:'= JJ rmO(^7^), (158) 

A=l,2,...,5 

where 

Tg'^ = exp(2 



B 







B 



32 



32 





(159a) 



r4'''' = exp(z 




^16 











032 



exp(iPo 



B 





is) 



-^t{^) 



B 



16 



(159b) 



16 
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r^'^ = exp(i 


" bI''' 






Bi'^ 




) = exp(zPo^ 






048 





B 



B 







(159c) 



= exp(i 


" 













056 



B 



B 



-^t(«) 







(159d) 



and 



r^''^ = exp(i 


" bI'-' 






B^^ 




) = exp(iP(f^ 






Oeo 





(159e) 



is defined by 

B^l = sinc(r^)e^4^ (160) 

for A e Zi,5. 

The above example generafizes as follows. For A = 1,2,3,..., s e -^0,2^-1-1 
and A e -^i,a, we have 

e^ti = J] r5f) + C»(y3^ , (161) 

A=1,2,...,A 

rSf) = exp(iPo®^-^ ® 

and 

S^l)=S(^sinc(:::^)e^^. (163) 
10.4 Expressing the ;52a 's in Bra-Ket Notation 

The CSD-rcady matrices T^^^ depend on matrices B^2>'- this section, we will express 

— (s) 

the B2X in bra-kct notation. 

Let us consider an example first, before dealing with the general case. Suppose 
A = 5 and consider the matrices Bgl^ that enter into the calculation of exp(i^g4 ). To 



(162) 
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calculate Sg^'' as defined by Eq. fll60p . we need to calculate sine 



^(0) 

calculate sinc(— |2-)e*~2^ using its Taylor expansion at x = 0, provided that we can 
calculate all powers of the matrix with reasonable efficiency. This can indeed 
be done because is a sparse matrix (it's effectively a 5 x 5 "tridiagonal band 
matrix"). Indeed, only rows and columns corresponding to nodes 1,3,7,15 and 31 are 
nonzero in .432^. It can be represented by: 



AO) 



(0) 

32 
2 



.(0) 



We can 



Hence, 



A 



(0) 
32 





1 


3 


7 


15 


31 


1 




di 








3 


di 




d2 






7 




d2 




ds 




15 






ds 




di 


31 








di 





(164) 



(-4S'; 



Other powers of ^32"* can be calculated just as easily. The pattern of entries that are 
zero in (.4.32'')" (for some integer n>l) can be easily understood by appealing to Dirac 
bra-ket notation. In that notation, it is clear that if a row number of (^32^)" is the 
starting point of a staggering drunk, then the column numbers of the non-vanishing 
entries of that row are the possible ffnal positions of the drunk. This, assuming that 
the drunk takes exactly n steps, either backwards or forwards, on the 5 state runway 
1, 3, 7, 15, 31. These rules imply that, at all times, the drunk will be (n minus an even 
number) steps away from his starting point. 





1 


3 


7 


15 


31 


1 







did2 






3 




d{ + rf^ 




(^2^3 




7 


did2 




dl + dl 




d^di 


15 




(^2^3 




dl + dl 




31 






d^di 







(165) 



Note that to calculate Sgg"*, we don't even need to calculate all the entries of 
32 ! (-^32^)" appears pre-multiplied by S32 



(0) 



we only need to calculate the last row of {A 



{0)^ 
32 y 



5^°^ - 4|31)(31|. This means that 
For example. 



-^32 (^^T = ^^5|31)(31|{ |31)rf4(15|15)4(7| + 131)^4(15115)4(311 }. (166) 

Thus, the drunk starts at 31. After 2 steps, he must be either: back at 31, or two 
steps away, at 7. 

The strand illustrations Figsl8|,[9|, and [7] can be used, in conjunction with our 
analogy to a drunken walker, to predict the general form of all that contribute 
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to exp(i^g4 ). One concludes that there must exisl0 some real numbers bj such that 
the following is true. For s = 0, 

^3? = |31)(&i(31| + .6i(15| + bf{7\ + zbf{3\ + 6f (1|) 

< = \l5)ml5\+^b\'{7\+bf{3\+^b\'{l\) 

<^ = |7)(&?(7|+^6^(3|+6?(1|) . (167) 

< = |3)(&|(3|+ .6^(11) 

= mm 



For s 



= m{bll{30\+^bll{U\+bf{6\+^bf{2\) 
Bf^ = \U){b\iaM+^bl<^m+bl'{2\) 

B^s'' = mm + ^bl{2\) ■ (168) 

B^^ = |2)(6^(2|) 



B^^ = 



For k G Zq^i (and s = 2 + k E ^2,3), 



B^^^'^ = |29 - k)ibll_,{29 -k\+ ^bf,_,{13 - k\ + bf_,{5 - k\) 
B^;'^ = |13 - k){b\U{13 -k\+ ^^,(5 - k\) 
Bt'^ = \5-k){bl,{5-k\) 



For k G Zo,3 (and s = 4 + G ^4^7) 



For A; G Zqj (and s = 8 + k E Z; 



8,15), 



B^^^^ = \23 - k){b^_f,{23 - k\) 

^16 ~ ^8 ~ ^4 ~ ^2 ~ 



(169) 



SS^'^ = |27 - fc)(6P_,(27 -k\+ blUill - k\) 

BT^ = \n-k){b\U{n-k\) . (170) 

_(4^.) ^ _(4+.) ^ -(4+.) ^ Q 



(171) 



The above example generalizes as follows. For A = 1, 2, 3, . . ., s G ^0,2^-1-1 
and A G Zi /^, we have 



^ An important special case is when the d\ (defined by Eq. (|13ip ) are all equal to V2g. Appendix 
IB] gives an explicit expression for the bj in this case. 
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b['^ ^ 12^ _ s _ 1) J2 ''""'K'^s-ii^'' - « - 1| ■ (172) 

Note that B^x = if A < A™„(s). 
10.5 Finding SVD of B^'^'s 

In this section, we find an SVD for each ■ Luckily, such SVD's are one sided, and 
can be found in closed form. 

As an example, assume A = 5, and consider all that enter into the calcu- 
lation of exp(i^64). 

Lemma 4 The following SVD holds for some pH, 6^1, Ofj, Of, Of e M; 
proof: 

In this proof, the symbols 6, p, should all have a superscript of 32. We will 
omit this superscript to simplify the notation. 
In the previous section, we showed that 



4? = |31)(63i(31| + ^6i5(15| + 67(71 + + 6i(l|) . (174) 



Note that 



31 = (11111), 15 = (01111), 7 =(00111), 3 =(00011), 1 = (00001) . (175) 
If we define pz and 0^ by 

Ps^Jbl + bl, cos^3 = -, sine3 = ^, (176) 
^ Ps Pa 



then 

{ih,bi) = (ip3,0) 

Thus, 



63 -ibi 
-ibi 63 



^ (zp3,0)e^^3^- . (177) 



Pa 



631(311+^615(151+67(71+ \_f 631(311 +i6i5(15| +67(71+ \ i^(4)n(a)n(2)e3ax(l) 

+i63(3| + 6i(l| ; \+^P3{^\ ) 

(178) 
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If we define pi and ^7 by 



p7 = + p'i , cos 6*7 = — , sin 6*7 = — , 



P7 



P7 



then 



(^7, ips) = (P7, 0) 



67 ips 
ip3 h 



-^{p,,0)e''^-- 

P7 



Thus, 



(179) 



(180) 



031 (31| +2615(151 + 
+br{7\+zps{3\ 



&3l(31| +i6i5(15|+ \ m(4)n(3)6'7<Tx(2) 
+P7(7| ' 



(181) 



If we define pis and 6*15 by 



2,2 /I ^15 • /) ~P7 



Pl5 = V + P7 > COS 6'i5 



Pl5 



Pl5 



then 



(182) 



{ibi5,p7) = {ipi5,0) 



bi5 -ip7 
-ip7 bi5 



— = (ipi5,0)e^^«'^^ . 

Pl5 



Thus, 



&3l(31| + 

+2615(151 + P7(7| 
If we define psi and 6*31 by 



631(311 + 
+^Pi5(15| 



^m(4)6liBCTx(3) 



(183) 



(184) 



/1.2 , 2 n "31 . /) Pl5 

P31 = V ^31 + Pi5 > COS 6'3i = — , sm 6*31 = — , 



P31 



P31 



then 



(&31,«Pl5) = (P31,0) 



631 ipib 
ipi5 hi 



P31 



(P3i,0)e 



Thus, 



(185) 



(186) 



( 631(311 +^Pl5(15| ) =p3i(31|e''3i'^-(^) . 



(187) 



QED 

The techniques used in the proof of the above claim can also be used to find 
an SVD of for all s G ^o,i5 ^ ^ -^1,5 (These are all the that enter into 
the calculation of exp(i^64))- One finds the following. 
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For s = 0, 



= P3l|3l) (31|e^^3?^x(4)gm(4)032<,^(3)gm(4)n(3)ef ax(2)gm(4)n(3)n^ 
Bjg = pl6|;^5^^;^5|gielfax(3)gm(3)0pax{2)gm(3)n(2)0l6^x(l) 

=p8|7)(7|e^^?'^^(2)gm(2)eiax{l) . (188) 

Bf = p||3)(3|e^^t-x{i) 



For s = 1, 



^32 = p^g|3O)(3O|e^^3O<^x(4)gm(4)032ax(3)gm(4)n(3)ef ax(2) 

Me = pil|14)(14|e^^i4'^^(3)gm(3)0i6ax(2) 

_p8|6^^6|gie|ax(2) . (189) 



=p^|2)(2| 

B^^^ = 



For k e Zo,i (and s = 2 + A; e ^2,3), 

Mr'^ = Pi-fc|29 - k){29 - k]e'(>'iU<^xi4)^mmi_,ax{3) 

Bf^'^=pU\5-k){5-k\ 



For A; e Zo,3 (and s = 4 + A; e ^4,7), 



Hir^=plt.|ll-A:)(ll-fc| ■ (191) 



For A; e Zo,7 (and s = 8 + A; e ^8,15)) 



-^32^*^^ =Pi-ifc|23-A;)(23-A;| 

^(8+.) ^ ^(8+.) ^ ^(8+.) ^ ^(8+.) ^ Q 



(190) 



(192) 



10.6 Expressing the T^^^'s as EG Circuits 

In this section, we express each F^*^ as an EG (elementary gates) circuit. 

Let us consider an example first, before dealing with the general case. Suppose 

A = 5. 

Consider s = 0. 
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Plugging the SVD of given by Eg. (11 88 p . into the definition of Fg , given 
by Eq. fll59ap . one finds 



ul 



32 



'32 



^iplj{\63){3l\+h.c.) 



u- 



32 



'32 



(193a) 



= [/32(4, 3, 2, 1, 0)t«(5)e^Pi?-^(5)"(4)'^(3)'^(2)r^(l)n(0)f;^^^4^ 3^ 2, 1, O)'^^^) ,(193b) 

where U32 is defined by 



f/32(4,3,2,l,0) = e*'^3i-xWg 
Eq. fll93bp simphfies to 



(194) 



^ g-m(5)n(4)n(3)n(2)0|Vx(l)g-m(5)n(4)n(3)0f ax(2)g-m^^^ 

0{e*'^3?'^x(5)n(4)n(3)n(2)n(l)n(O) j " 

(195) 

Plugging the SVD of i^lg, given by Eq. (I188P , into the definition of r4 , given 
by Eq. (]159bp . one finds 



.(0) 



U- 



U- 



-\ 02 



16 



'16 



(196a) 



16 „iPo®m|31>(15|+?i.c.) 
/l6j 

= Uie{3, 2, 1, 0)t"(4)eM5)p;i-x(4)n(3)n(2)n(l)„(0)f;^g^3^ 2, 1, 0)"^^^ , (196b) 

where Uiq is defined by 

;7l6(3, 2, 1, 0) = e*^l5^-^(3)gm(3)elVx{2)gin(3)n(2)el6ax(l) _ (;L97) 

Eq. (1196bp simplifies to 



(198) 

Similarly, one finds 



r 



(0) 



^iP®2®(p8|15>(7|+h.c.) 



(199a) 



Usi2, 1, 0)t"(3)eM5)"(4)pf-x(3)n(2)n(l)n(0)fjg^2, 1, 0)"^=^) , (199b) 



where 
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[/8(2,1,0) = e^^?'^^^'^e 



Eq. (fT99b|) simplifies to 



T(3)n(2)0|(7x{l)g-m(3)e»<7x(2) q |gm(5)n(4)p8<7x (3)n(2)n(l)n(0) | 



Similarly, one finds 



(200) 



(201) 



.(0) 



,iP®2^(pf|15>(7|+h.c.) 



where 



t/4(l,0) = e''3-x(i) _ 

Eq. fl202b|) simplifies to 

Similarly, one finds 



(0) 



,iPo®4®(p?|3>(l|+ft.c.) 
^m(5)n(4)r!;(3)n(2)pf<Tx (l)n(O) 



To summarize, we have shown that: 



(202a) 
(202b) 

(203) 
(204) 



(205a) 
(205b) 



.(0) 
4 

^(0) 



-m{5)n(4)n(3)W(2)9f ax{l)p-m(5)n(4)W(3)9Pax(2)„-m^ 



ax (5)n(4)n(3)n(2)n(l)n(0) 



} 



.(0) 



-m(4)n(3)n(2)9i6^x(l)g-m{4)n(3)9l6^x{2)g-m(4)elf q |gm(5)pJ6f7x(4)n(3)n(2)n(l)n(0) | 

-m(3)n(2)e|(Tx{l)g-m(3)6l»crx(2) q |gm(5)n(4)p»(7x (3)n(2)n(l)n(0) j 
-m(2)0|c.x(l) |gm(5)n(4)n(3)p|f7x(2)n(l)n(0) j 



(0) _ gin(5)ri:(4)n(3)ri:(2)p2crx(l)n(0) 



(206) 



FigHni represents the SEOs of Eq. fl206p as circuits. 

Using the same method that we used for s = 0, we find for s = 1: 



r^^^ = e-*"(5)n(4)n(3)9|Vx(2)g-m{5)n(4)932^^(3)^_i7T(5)eM ^ jgiplg^x (5)n(4)n(3)n(2)n(l)n(0) | 

^ g-m(4)n(3)0f <7x(2)g-m(4)e;|(7x(3) q |gjn(5)p;|fTx (4)n(3)n(2)n(l)n(0) j 

pW ^ g-m(3)e|crx(2) |gm(5)n(4)p|(7x(3)n(2)n(l)n(0) j 

_ gm(5)n(4)n(3)p|o-x(2)n(l)n(0) 

r« = 1 

(207) 
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Figure 10: T^'^'s for s = 0. 



FigHH represents the SEOs of Eq. (120 7p as circuits. 
For s = 2, we find 

p(2) ^ g-m(5)n(4)ef|<7x{3)g-m(5)eicrx{4) q j-gip^l^x (5)n(4)n(3)n(2)n(l)n(0) | 
T-(2) _ -m(4)9lfax(3) r m(5)p}|f7x(4)n(3)n(2)n(l)n(0)| 

i -* ('?C\R'\ 

p(2) ^ gjn(5)n(4)p|<7x(3)n(2)n(l)n(0) ' 

rf = ) = 1 

The rlj*"* for s = 3 are identical to those for s = 2, except that one must replace n(0) 
by rT(0) in them. The angle parameters (i.e., pi,Oi,^k) ^^^o change. FigJT2] represent 
the SEOs of Eg. (12081) as circuits. 
For s = 4, we find 

p^4) ^ g-m(5)032^^(4) Q |gip32^^(5)„(4)„(3)^(2)„(i)„(o)| 

r^4) ^ giW(4)p;fax(4)n(3)n(2)n(l)n(0) _ (209) 

r(^) = r(^) = ) = i 

The r^^^ for s = 5, 6, 7 are identical to those for s = 4, except that we must re- 
place n{l)n{0) by n(l)n(0), n(l)n(0), and n(l)n(0), respectively, in them. The angle 
parameters also change. FigHni represents the SEOs of Eq. (12091) as circuits. 
For s = 8, we find 



f(8) _ g«P23<^^(5)"(4)n(3)n(2)n(l)n(0) 



r(8) _ p(8) _ p(8) _ p(8) _ . • (210) 

■'-4 — ■'-3 — 2 — 1 — ^ 

The rl^'' for s = 9, 10, . . . , 15 are identical to those for s = 8, except that we must 
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r,(i) 



Figure 11: T^'^^'s for s = 1. r\'^' is not shown because it equals 1. 



^(1) 



replace n{2)n{l)n{0) by the appropriate projection operator in them. The angle 
parameters also change. Fig{Tl] represents the SEOs of Eg. (12101) as circuits. 

The above example generalizes as follows. Suppose A = 1,2,3,... and s G 
■^0,2^-1-1- Then T^*^ = if A < Xmin{s) since for such A, = 0. For A G Za„„,(s),a, 
there exist (p^^^j , 6^^^ G M such that 



^ J-j- |gm(A)n(A-l)...n(A-j>(^^^,ax(A-j-l)| 

j=A-A„i„(s)-l, ...,2,1,0 

Q|gm(A)n(A-l)...n(A+l)eWax(A)P5-^(A-l,...,l,0)| _ ^211) 

In Eq. (l21ip . we use bin{s) = (sa_iSa-2 • • • siSq) to denote the binary representation 
of s, and bin{s) = (sa-iSa-2 • • -SiSq) to denote its bitwise negation. 



10.7 Weaving the Strands 

In this section, we finally achieve our goal of compiling e*^"". 

Suppose a and /? are n-tuples (ordered sets), not necessarily of the same length, 
of distinct qubit positions. Assume that these two n-tuples are disjoint (if treated as 
non-ordered sets). Let vr and tt' be commuting projection operators acting on qubits 
p. Let M be any operator acting on qubits a. A simple identity that will be useful 
in what follows is 

M{ay^f^^M{af^'^^ = M {aY'^^^'^^' . (212) 

This identity allows us to combine operations that are being needlessly performed 
separately. 
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r(2) 

^3 



-Q- 



-6-Rx-6- 



o- 



r(2) 



-o- 



r,(2) 



-o- 

-Q- 
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-o- 



-o- 



-o- 



-o- 



-o- 



-rTU-Ir" 



-cHR^n-6-'-6|R 



o- 




r,(3) 



Figure 12: r^*'''s for s = 2, 3. T^^ and Fg*'' are not shown because they equal 1. 



Define F^'^^ for s G -^0,2^-1-1 by 

Yis) = emVi = Yl T^^^ + 0{g^) , (213) 

A=1,2,...,A 

where we have used Eg. (11611) . 

Note that the F*-**^ for different strands commute, because any two different 
strands "five" on disjoint sets of nodes. The main point of this section is that many 
operations can be combined if we multiply the F'^'^^'s judiciously. 

For definiteness, consider our usual example of A = 5. 

Note from FiglHl that strands can be grouped into equivalence classes {0}, {1}, 
^2,3, and .^8,15- AH members of the same equivalence class start and end at the 
same tree levels. They also have the same strengths diS for transitions between levels. 
One could say that they are identical in every respect except that the names of their 
nodes differ. 

Now consider Fig. fini) . which gives circuits for the F^^^'s, for strand class ^8,i5- 
This figure does not tell us what are the specific angles for the qubit rotations of 
the type R^. But from our observation that equivalent strands are "identical except 
for node relabelling" , it is clear that any qubit rotation of type Rx in a strand si is 
identical to the "analogous" one in any strand S2 of the same class. Therefore, we 
can apply identity Eq. (l212p to combine the F*^'^^ of all the strands of class ^8,15 • 
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Figure 13: Fj^ s for s = 4 to 7. F]^ , Fj and Fg are not shown because they equal 
1. 



n r 



(214a) 





Likewise, by applying identity Eq. fl212p to Fig. (fT3l) . one finds that 



(^) 



(214b) 



Finally, by applying identity Eq. (12121) to Fig. lIT^ . one finds 
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Figure 14: r^*'''s for s = 8 to 15. Tf\ ^t'^ shown because they 

equal 1. 



-Rx 



-Rx 
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. (214c) 



10.8 Complexity of e*^*'^ Compilation 

In this section, we will show that our SEO for e^^^'' contains (9(A^) CNOTs, in the 
case of balanced binary NAND trees. 

If S represents a SEO or its associated circuit, let C{S) denote its circuit 
complexity. We will evaluate the circuit complexity of our SEO for e*^*"^ by counting 
the number of control vertices in its quantum circuit. 

Assume a tree with 2'^"'"^ nodes. According to Eq. fll40p . 



C{e'^'^) = 2C(J'2A+i) +C(exp(U2A+i)) • 
According to Eq. (11381) . 

A 

C(^2A+i ) = ^ C(F2. (A - 1, . . . , 2, 1, 0)"W) . 

A=l 



(215) 



(216) 
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An exchange operator E4{a, [3) can be expressed as product of 3 CNOTs. So a singly- 
controlled exchange operator has six control vertices. This observation and Eg. (11301) 
imply that 



C(F2"/^)) = 6A-5. (217) 
Combining Eqs.([2l6]) and (I2T7D yields 

C(J-2A+i) = C(A') . (218) 
According to Eqs.([Iin]) and fl^ . 



c(exp(zAA+0) = c( n r(^^) = c(r(°))+ c( J] • (219) 

s=0 o-eZo,A_2 s=2<^,2<^+l,...2<'+i-l 

Here a labels the different strand equivalence classes. According to Eqs. (l214p . for 

O" G ^o,A-2, 

c( n rW) = E{C(rr)-^}. (220) 

s=2<^,2<^+l,...2''+i-l A=l 

According to Eq.([2lID and FigsHO] to [H 

C(rl^)) = A + 2 5^ {j + 1}. (221) 

j=0 

Combining Eqs. (I219p . (I220p . (I22ip . and doing some algebra, yields 

C(exp(M2A+i)) = C(A^) . (222) 
Now we can use Eq.([2I5D, fl^THD and fl22^ . to conclude that 

C(e*^*'-) = O(A^) . (223) 



11 Compiling Input Graph ("Oracle") 

In this section, we will show how to compile the evolution operator e*^™ for an input 
graph ("oracle"). 
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Suppose A = 1, 2, 3 . . .. Consider a tree with Ni^g = 2 leaves, with leaf inputs 
Xk G Bool for k G Zq 2^_i. The input Hamiltonian for such a tree is: 



Hi, 





Xo 

Xi 




Xo 








Xi 

















where G M, and is defined to equal x^^^^j^y 
Define operators U and vr by 

Let 



(224a) 



(224b) 



(225) 



e = (A-l,...,l,0). 



(226) 



Then 



^iHi^ = [[/(A)]-(«) (227a) 
= ax(a)"(«^)f/(A)"(")(Tx(a)"^«^)|0), . (227b) 

Here a is an ancilla qubit. To get Eq. fl227bp . we did the same as we did in Eq. (l40|) : 
we replaced the multiply controlled U (2) operator by a singly controlled one, and two 
MCNOTs. So, to compile e*^'", all we need to do is to compile ax{ciY'^^'^ ■ 

The operator ax{ciY'^^^ appears frequently in the quantum computation liter- 
ature, but usually in a different guise. It's more common to denote the function x, 
in Eq. (12251) by a more standard name for a function, such as /a : Bool^ — > Bool. In 
this notation, the operator vr of Eq. (12251) becomes 

vr= 5^ /a(x')P-, . (228) 

x'(^Bool'^ 

Thus, 
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ax{ap^^y)a\x)^-= \y®fA{x))^\x)^-. (229) 

The best general way of compiling axiaf^^ is not known. It is known that there are 
some /a() for which crx(a)'^*-^^ has complexity 0(2^) = 0{Niys). In Appendix O we 
give a compilation for a special case, that of a "banded oracle" (an oracle for which 
the ordered set of NAND formula inputs (xfe)vfe has a fixed number of bands of ones). 
Our compilation of the banded oracle evolution operator has complexity (9(A^). 

12 Overall Circuit Complexity 

The following table summarizes the circuit complexity of the various compilations 
that have been presented in this paper: 



evolution 
operator 


description 


number of CNOTs in SEO 


exact? 




balanced binary NAND tree 
with 2^ = A'";^,^ leaves 




no 




loop of length Ns 


0{Ns) 

0{NiJ) for Ns ^ NiJ 


yes 




glue connecting tree and loop 


OiA) 


yes 




cuts a loop open 


0(A) 


yes 




line of length Ns 


0{Nb'') 


no 




oracle, with 2^ = Ni^s inputs 


Could reach 0(A'/^,s). 
(9(A^) for banded oracle. 





Henceforth, we will say that a SEO is tractable if it has 0(A'^) CNOTs, for 
some k > 0. 

Our plan is to combine e'^^''^"^ and e'^^""'"'' via Trotterized Suzuki. Let Nexp be 
the number of times e*^'"'"= or e^^""'"'' appear in the Trotterized Suzuki product that 
approximates e^^^'"^"'~^^''°''''\ According to Eg. (15^ . if we want to evaluate a NAND 
formula in time ~ 0{{Niys)^~^'^), then 

N,^p^O{{Ni,s)'^+') . (230) 

Thus 

(^(^^i{H,^ik+H,orr)-^-^ ^ A'e^pC(e^-^'"'''=e'-^'=°") (231a) 

^ 0{{Ni^s)^+')C{e'^'"^''e'^^''"-) . (231b) 

If f,iHtuik(,iHcorr ig tractable, then C(e*(-^'"'''=+^"°'-'-)) = 0{{Ni^s)^+^)- 



54 



A Appendix: Exponential of Anti-block Diagonal 
Matrix 



Lemma 5 Suppose F 


e C 


ixn 


Then 




' 


Ft " 




exp(i 


F 





) = 



cos 



(\/FtF) 



■ sin(VFFt) 



„. sin(VFtF) pt 



COS 



(232) 



proof: 



Note that 



Since a+ = |0)(1| and cr_ 
hand side of Eq. ( 



Ft ■ 
F 

= |i)(o|, 



(g) F + a, (g) FK 



(233) 



= (T-^ = 0. Thus, when we raise the right 
to a power, only terms proportional to an alternating sequence 
of and cr_ survive. For n = 1, 2, 3, . . ., one gets 



and 



(a_ ® F + a+ ® Ft)2" = (a+a_)"® (FtF)" + (a_a+)"® (FFt)" (234a) 

= Po® {F^Ff + Pi® {FF^Y , (234b) 



(a_ ® F + (T+ ® Ft)2"+i = Paa+®{F^FYF^ + Pia^®{FF^YF (235a) 

= a+ ® (FtF)"Ft + a_ ® (FFt)"F . (235b) 



Thus 



n=0 ^ 



2n 



= Po ® cos(\/FtF) + Pi ® cos(yFFt) , (236) 

and 



E 

n=0 



(2n + 1)! 



[a_ ® F + a+ ® Ft) 



2n+l 



0-4 



sin(VFtF) 
z \ ' Ft _^ ^_ 



VfTf 



.sin(VFFt) 
i , — F. 



VFFt 



(237) 



QED 
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B Appendix: Special Case dx = \/2g for all A 

An important special case is when dx = V^g for all A, where dx is defined by Eq.f ll3ip . 
This appendix will discuss some idiosyncrasies of this special case. 

When the dx are all equal, it is possible to make a slight redefinition of in 
Eq. f illip so that the decomposition of e^'^a'^+i into given by Eq. flllOp is good to 
order 0{g^) instead of 0{g^). We show this next for A = 3. Generalization to higher 
A will be obvious. 



Lemma 6 Let 



Then 



where 



/ = 1 + ^ 



gMi6 = r<3>r<3>r<3> ^ ^^g4 



^<3> 



exp(« 



_0 Bs 
B, 












B^ 




P<3> 
J- 2 


= exp 




B, 

















Os J 











B2 




p<3> 
J- 1 


= exp( 




B2 















0l2 J 



(238) 

(239) 

(240a) 
(240b) 

(240c) 



and 



B2>- = B2X sinc(^^^)e* 2^ , for A G ^1,3 



(241) 



proof: According to Eq. (l67|) . 



exp(« 



exp(? 



^8 ^8 

Bs 

Or 



exp( 



T1 + T2 
^Bs 



iB. 



—2 



(242) 



0(/) . (243) 



One finds that 
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^4 



T2 = 0{g' 



(244) 



and 



BsAsBs = 



(245) 



Therefore, 



JA16 



exp[i 



As 



0« 



exp[i 



^4 





i38 

08 



fB, 



)Tf> + 0{g') 



(246) 
(247) 



To go from Eq. fl2i6|l to Eq. fl27n) . we used A 



A -64 




Furthermore, 



exp(z 



/S4 



exp(z 



Aa 



exp(z 



6 



0_ fB2 
fB2_ 

/i34 



/S4 
O4 



P<3>p<3> 
■■-1 '-2 ) 



where we used A4 
QED 



1248) 
(249) 



A2 
B2 



B2 




and g'^f ^ g"^ 



Another useful consequence of dx = V^g for all A is that in this case it is easy 
to calculate the 6^ coefficients used in Section 110.41 Indeed, observe that in this case 
all non-zero entries of equal V^g- Define 



g\/2 smc(^-^)e'~ 



(250) 



Let hr and hi denote the real and imaginary parts of 6, respectively. Now we can use 



the following rule. In Section [10.41 in the bra-ket expansions of Sgl^ replace, all 6^ 
that are preceded by an imaginary i, by hi , and replace, all h^ that are NOT preceded 
by an imaginary i, by hr- For example, in the example of Lemma HJ one has 
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P3 = = l&l , cos^3=|^, sin^3 = 7^, (251) 



P7 = 


Vb'r + \b\' 


, COS 6*7 = — , 
P7 


sin 6*7 = 


P3 
P7 ' 


(252) 


Pl5 = 1 


^bj + P? , 


COS 015 = , 

Pl5 


sin^^is = 


-P7 
Pl5 


(253) 


P31 = ' 


\lbl + pl^ , 


. br 
COSt^si = — 

P31 


, sm6'3i = 


_ Pl5 
P31 ' 


(254) 



C Appendix: Compiling Banded Oracle 

Suppose G Bool for b G Bool^, and A^^/^js = 2^. Suppose a is a qubit, and 
,^ = (A — 1, . . . , 1, 0) is a vector of A qubits, and all these A + 1 qubits are distinct. 
In this Appendix, we will show how to compile the oracle evolution operator 

Uo = (Tx{a)^^b^BooiA^sPs(i) ^ (255) 

when X = [xq, xi, . . . , xni^^-i]'^ is banded. If we envision x as a sequence of ones and 
zeros, then we will call a band of x, any subsequence of x consisting of adjacent terms, 
all of which are one. It x contains a fixed (A independent) number of bands, we will 
say it is banded. 

We will consider first the case when x has a single front band. By this we 
mean that the first A^^i terms of x are one, and the rest are zero. Thus 

X = ^^1^^^, 0^0^^^ f ■ (256) 

TVi times Niyg—Ni times 

Define a binary vector bmax and its corresponding set of binary vectors S(bmax) by 



bmax = bin{Ni - 1) , S(bmax) = {b E Bool^ : < dec(h) < decipmax)} ■ (257) 
In the single front band case, Uq reduces to: 

Uo = (Txia)^'^^'-'----''^'^^'^ . (258) 

Uo, in the form given by Eq. fl258p . is already compiled. But the length of this 
compilation can be reduced significantly by reducing J2beS{b„T,ax) ^b- 
sider an example first, before dealing with the general case. Suppose A = 8 and 
bmax = (0110, 1101) = 6m(109). Although Zl6es(fem(io9)) is a sum of 109 projection 
operators of the type where b G Bool^, it can be expressed as a sum of just six 
simpler projection operators: 
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fee5(fem(109)) 



Poo--,---- 

+-Poio-,---- 

+-Poiio,o--- 

+-Poiio,io-- 

+-Poiio,iioo 

+-Poiio,iioi 



(sum from 0000,0000 to 0011,1111) 
(sum from 0100,0000 to 0101,1111) 
(sum from 0110,0000 to 0110,0111) 
(sum from 0110,1000 to 0110,1011) 



(259) 



FigHSl gives a circuit diagram for Eg. (12581) . assuming dec{bmax) = 109, and with 
X]b65(fem(io9)) expressed in the simphfied form given by the right hand side of 
Eq.( 



Uo = 







-1 

-0 1 

-1 2 

-1 3 

-0- 

-1- 



(*) 



O — (I — It — ly 



O — (I — II — () — 1^ 



Figure 15: Circuit diagram for Eq. fl258p . when the first 109 components of x equal 
one, and the final 2* — 109 = 256 — 109 = 147 components equal zero. 



It's also interesting to consider an example in which dec{bmax) is an even 
number instead of an odd one. When decipmax) = 108 instead of 109, we must 
remove the operator marked by an asterisk in Eq. fl259p . and the operator marked by 
an asterisk in FigJTSl 

The pattern of the control vertices in Eq. (12591) and FigHSl is not hard to 
uncover. There is exactly one MCNOT for each nonzero bit in 6m(109). In addition, 
there is one "final" MCNOT with controls equal to Pbin{im)- Call 1C{(3) the "non- 
final" MCNOTs. There is precisely one of these for each bit (3 G Zq^\ such that the 
(5 component of hmax (i-e., (^max)/?) equals one, and none when {hmax)f3 = 0. /C(/5) 
always has control Pq at bit jS. At bits a G Z/3+i^a_i, /C(/3) has a control P(%^^^)^- 
Hence, in general. 
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Note that the number of MCNOTs in Uq is C(A). The number of CNOTs in 
each of these MCNOTs is C(A). Thus, the number of CNOTs in Uq is ^(A^). 

Now suppose X has a single band which is, however, not at the front. Suppose 
it ranges from dec{bmin) 7^ to decipmax)- Then just multiply Uq for a single front 
band up to dec{hmin) — 1 times Uq for a single front band up to dec{bmax)- Multiple 
bands can be handled similarly. 

Sometimes it is possible to apply ax on individual qubits and/or apply qubit 
permutations to Uq so as to get a new Uq with fewer bands. Fewer bands will lead to 
a shorter SEO of the type proposed in this appendix, ax on individual qubits and/or 
qubit permutation^ do not increase the length of a SEO if they are applied to the 
SEO on both sides, via a product. 



D Appendix: Compiling Evolution Operators in 
Grover's Algorithm 

The goal of this appendix is not to say something new about Grover's algorithm [27] . 
After all, Grover's algorithm has been studied so extensively in the literature that it's 
almost impossible to say anything new about it. The goal of this appendix is, rather, 
to review how one compiles the oracle and non-oracle evolution operators associated 
with Grover's algorithm. This will allow the reader to compare the compilation of 
Grover's algorithm with the compilation of FGG07 presented in this paper. 

Let Nb be the number of qubits and A^^ = 2^-^ the number of states used in 
Grover's algorithm. Let Xo € Bool^^ be the target state of the algorithm. Define 



1 



xeBool^B 



(261) 



and 



1 1 

1 1 

1 1 

1 1 



(262) 



Thus, in matrix notation, \fi) is an Ns dimensional column vector, and fi is an Ns x Ns 
matrix. 



^Don't confuse qubit permutations with state permutations. 
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In Grovcr's algorithm, one alternates between an oracle evolution operator, 
call it Ucorr, and a non-oracle evolution operator, call it Ubuik- The oracle evolution 
operator depends on the target state whereas the non-oracle evolution operator does 
not. More precisely, these two evolution operators are defined by: 

U^rr = = (-1)^^'' = e'""^^" , (263) 

and 

Uhuik = 2|/x)(/x| - 1 = = e*-(^+i) . (264) 

In Grover's algorithm, one applies ^/Ns times the product UcorrUbuik- A variant 

1 1 

of this would be to apply Nt = Ns^~^^ times the product Uc^rr U,^fi^ , for some 6 > 0. 
One can think of this variant of Grover's algorithm as a Trotterized Lie approximation: 

AT" 

e J^s^+'e j , (265) 

for some f e R, where the Hamiltonians Hcorr and H^uik are given by 

(266) 





_ 7-/^ _ 
— <^corr — 


Hcorr 




^bulk 




Ns" 


1 








_ TT^ _ 




_ 7r(l+M) 


~ '^bulk ~ 




Ns" 



This variant of Grover's algorithm is a quantum walk over a fully connected graph 
with Ns nodes. Transitions occur on this graph along edges connecting distinct nodes 
and also from a node back to itself. Self transitions occur with strength proportional 
to 1 -I- for all nodes except the one representing the target state. Self transitions for 
the target node occur with strength proportional to 2 -|- larger than the strength 
for the other self transitions. Transitions along the edges connecting distinct nodes 
occur with vanishing strength proportional to 

— * 

Compiling the oracle evolution operator for Grover's algorithm is trivial. Let ^ 
denote the iVs-dimensional vector of qubit positions for the Nb primary qubits used 
in Grover's algorithm, and let a denote the qubit position of an additional (not in ^) 
ancilla qubit. Then 



Ucorr = {-If '"^^^ = (-l)"(")^-«"^|l)a (267a) 

= a^(a)^-«)|l)« (267b) 

- H{a)ax{af^''^^\Hax){a)\0)a . (267c) 

We used az = (-1)", (7x|0) = |1) and HaxH = az- 
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Compiling the non-oracle evolution operator for Grover's algorithm is also 
trivial. One notes that 

I;,) =i7«^i3|o) ^ (268) 
since the first column of H®^'^ is all ones. Therefore, 

^_^®iVs |o)(o| . (269) 

Eg. fl269p is merely the eigenvalue decomposition (and SVD) of /i. Since /i is a circu- 
lant matrix, we could have obtained Eg. fl269p from the eigenvalue decomposition of 
circulant matrices presented in Section [STTl An immediate conseguence of Eg. (12691) 
is 



U,^ik = (_i)lo)(o| H®^^ = -H^""^ (-l)^«-° « if^^^ . (270) 

One can compile (— l)-'^o.o. ■•.o in the same way that (— 1)'^*° was compiled above. 

Note that we have expressed both Ubuik and Ucorr as SEOs containing a single 
MCNOT. This MCNOT can be expressed as a SEO with 0{Nb) CNOTs. 
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